Sub-gap conductance in ferromagnetic-superconducting mesoscopic structures 
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We study the sub-gap conductance of a ferromagnetic mesoscopic region attached to a ferromag- 
netic and a superconducting electrode by means of tunnel junctions. In the absence of the exchange 
field, the ratio r = 7/eT of the two tunnel junction resistances determines the behaviour of the 
sub-gap conductance which possesses a zero-bias peak for r ^> 1 and for r <C 1 a peak at finite 
voltage. We show that the inclusion of the exchange field leads to a peak splitting for r < 1, while 
it shifts the zero-bias anomaly to finite voltages for r> 1. 

PACS numbers: 74.50 +r, 72.15 Nj 
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In the last few years there has been a great deal of 
activity in the study of transport through hybrid super- 
conducting mesoscopic structures (S/N) ]IJ. Amongst 
the many interesting findings it is worth mentioning the 
observation of a sub-gap conductance in S/N struc- 
tures, long-range || and anomalous proximity effects ||] 
in N-S systems and giant oscillations in the electrical 
conductance ||^| and thermopower in S/N structures 
containing superconducting and/or normal loops. 

In parallel with these discoveries, ferromagnetic- 
normal hybrid systems (F/N) have stimulated the cu- 
riosity of the scientific community || and more recently 
combined structures involving both ferromagnetic and 
superconducting materials like F/S or F/N/S have been 
studied both theoretically PJIo| and experimentally (jll) . 
We note here that, in general, the introduction of the 
exchange field leads to a decrease of the Andreev scat- 
tering due to the change in the ratio of the two spin 
populations. For our case, where the exchange energy is 
of order the superconducting energy gap, the density of 
states is almost constant and this reduction effect will be 
small. 

In this paper we study how the presence of a finite ex- 
change field affects well-known coherent phenomena in 
disordered N(F)/S systems. We shall focus attention 
upon the F/S structure shown in Fig.l, calculating its 
electrical conductance in the diffusive regime by means 
of both analytical and numerical techniques. The former 
is based on the use of the quasi-classical Green's func- 
tion approach, whereas the latter solves the Bogoliubov 
de Gennes equations on a tight-binding lattice. Both of 
these methods have been used extensively in the study of 
N/S systems M. For our structure, we recall that in the 
ballistic limit for a simple NIS system the conductance is 
suppressed at voltages below the superconducting energy 
gap and the conductance-voltage curve is featureless. A 
sub-gap anomaly in the conductance appears only in the 
presence of two barriers. In the diffusive regime the sub- 
gap conductance exhibits a peak below the superconduct- 



ing energy gap A even in the absence of a second barrier, 
because the disorder itself provides the mechanism for 
multiple scattering. The position of this peak depends 
on the particular geometry and relevant parameters |L2| . 
In the absence of the exchange field, the ratio r = j/et 
of the two tunnel junction resistances controls different 
regimes for the sub-gap conductance: a zero-bias peak is 
present for r>l and a finite bias peak for r <C 1 . Here 7 
and ct are proportional to the conductances of the N'/N 
and N/S interfaces respectively. In what follows we show 
that the inclusion of the exchange field (we denote the 
exchange energy by h) leads to a peak splitting for r <C 1 
for small h (weak ferromagnetism) while the zero-bias 
anomaly is shifted to a finite voltage eV = h for r> 1. 

To numerically model the structure of Fig.|l| we con- 
sider a dirty normal or ferromagnetic region (denoted N 
or F) connected to external reservoirs by a clean semi- 
infinite normal or ferromagnetic lead (denoted N' or F') 
and a clean superconductor (denoted by S). For our nu- 
merical calculation we use a tight-binding version of the 
Bogoliubov-de Gennes Hamiltonian: 
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where a = \(— 1) for spin-t(|). In this equation Hq = 
J2i N) e «(*l ~ ^S(y) l*)0'l i s tne standard single-particle 
Anderson model with (ij) denoting pairs of nearest neigh- 
bour sites, h = Y^,i models the exchange energy 
and A = |i)Aj(z| is the superconducting order pa- 
rameter. The width of the whole structure is W and the 
length of the disordered ferromagnet is d (in units of the 
lattice constant a). The off-diagonal matrix elements t, 
which determine the width of the energy band, are equal 
to 1 throughout the system except at the F/S and F' /F 
interfaces where the value of t — tps(F'F) is chosen to 
model a barrier. Such barriers may be due to either the 
presence of an insulating layer or to a mismatch of the 
electronic parameters of the adjacent materials. Within 
the ferromagnetic lead and the superconducting region, 
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the diagonal matrix elements are = eo (with eo = 1, 
which keeps the Fermi level away from the van Hove sin- 
gularity in the band centre). For simulating bulk disor- 
der in the ferromagnet, the ej in the scattering region are 
taken uniformly from the interval — U/2 < e^ — eo < U/2 
where U is the disorder strength. A numerical decimation 
technique G3] is used to compute the Greens function for 
each realization of disorder. From this one obtains the 
Andreev reflection coefficients Ra aa and the conductance 
G = R$ + Rf @, in units of fi. 

Our analytical treatment is based on the transport 
theory described in fL2fl . There equations for the quasi- 
classical Green's functions were employed in the diffusive 
regime with the corresponding boundary conditions at 
the interfaces [|l5| . We begin by recalling the main results 
for non ferromagnetic structures. The normalised con- 
ductance S(V) = (dI/dV)/G n (where G n is the conduc- 
tance of the isolated normal wire) of the system N' /N/S 
(exchange energy h is zero) is equal to (see [O) 



S{V) 



deF' 



m 



- X (e) 



Vd{e) 



£t M(e)eT 



(2) 



Here F' v = (/3/2)[cosh~ 2 (e + eV)0 + cosh" 2 (e - eV)0\ 
is the derivative of the distribution function, V is the 
voltage difference between reservoirs, (3 = 1/2T, i>d(e) = 
Re((e + + «7) 2 — e|0 is the density-of-states in 

the central TV region, 7 = pD/2RNad, et = pD/2Rsnd, 
p and D are the specific resistivity and diffusion constant 
in the middle region, d is the length of the middle region, 
Rn,sh are the N'/N or N/S resistances per unit area in 
the normal state and €d — D/d 2 is the Thouless energy. 
We assume here that the length d is small enough and 
the condition » T is fulfilled (short contact). The co- 
efficient m(e) determines the conductance of the normal 
wire and is given by 
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The coefficient M(e) in eq.(|2|) determines the energy de- 
pendence of the N/S interface conductance. For exam- 
ple, from eq.(107) of Ref. jjj, we have 
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Here A and T are the gap and the damping rate in the 
superconductor, respectively. The first term in eq.(^) is 
due to the quasi-particle contribution (at zero tempera- 
ture it is not zero only if eV > A). The second term 
in M(e) gives the so-called interference current in SIS 
Josephson junctions and leads to a sub-gap conductance 
in SIN contacts. 



The three terms in square brackets in eq.(g) represent 
the normalized resistances: the first term is the resistance 
of the central N region, the second is the resistance of 
the N'/N interface and the third that of the N/S in- 
terface. The quantity er is a pseudo-gap induced in the 
middle region by the proximity effect. One can see from 
the expression for Vd that in the case of a small 7 (high 
resistance of the N'/N interface ) V4 has a form typi- 
cal for a superconductor with the energy gap ex- If 7 
is large compared with ex, the singularity in the density 
of the states in the N film disappears due to the strong 
damping caused by the N' electrode. In Fig. 2 we plot the 
dependence of the normalized conductance S(V) on the 
applied voltage V at zero temperature assuming that the 
resistance of the N region is small (e^/ey — 10). One can 
see that in the case of large 7, the main contribution to 
S(V) is due to the last term in brackets (N/S interface 
resistance) and the sub-gap conductance has a zero-bias 
peak. In this case there is a gapless state in the central 
N region. When 7 is decreased the main contribution to 
S(V) is given by the second term in brackets; the sub-gap 
appears in the central N region and the peak is shifted 
to ex- the conductance is just that expected for tunnel- 
ing into a superconductor with a energy gap ex- Similar 
behaviour of the sub-gap conductance was discussed in 
Ref. g§. 

The numerical results in the absence of an exchange 
field are shown in Fig.^|. To produce this curve we chose 
W = d = 20 sites. The barrier at the N/S interface 
was chosen by setting t^s = 0.4. The different curves 
correspond to ijv'Ar = 0.0, 0.1, . . . , 1.0. The result was 
obtained by averaging over 500 realisations of the ran- 
dom potential with U — 2 (elastic mean free path is 8.5). 
One can clearly identify the transition from finite to zero- 
bias peak. One should note that for these parameters we 
cannot reach the regime of 7 S> £t even in the absence 
of the barrier at the N'/N interface and so we do not 
obtain the Lorenzian line shape expected in that limit. 
This can easily be obtained, however, by using a smaller 
value for t^s- 

Consider now the case when the central region is a fer- 
romagnet. We describe the ferromagnet by adding to the 
Hamiltonian an exchange term H exc = h^2 k (c\, ^c^f — 

c\ \Ck,i)- One can write the quasi-classical equations for 
the Green's functions taking into account this term. One 
can easily show [0 that the energy should be replaced 
by (e — h) in all terms corresponding to the ferromag- 
netic regions. Therefore eq.(||) remains unchanged if we 
replace e by (e — h) in all such terms (the terms corre- 
sponding to the condensate functions in the supercon- 
ductor are unchanged as we do not consider an exchange 
field in these regions). Alternatively one can obtain the 
same expressions for the F'/F and F/S interface resis- 
tances as in eq.(||) using the standard tunnel Hamiltonian 
H T = U T J2k, q , a ( c k, a a q,a + a \, a c k,a), where a 9iQ are the 
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destruction operators of the reservoirs. We find the con- 
densate functions induced in the ferromagnet due 
to the proximity effect 
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where e T = N s \U F /s\ 2 and 7 = N F >\U J 
the density of states in the 5 and F 1 reservoirs in the 
normal state. Using eq.(|^) we come again to the expres- 
sions for the resistances of the F/S and F 1 / F interfaces 
(the third and second terms in eq.(^)). In Fig. [| we 
present the voltage dependence of the conductance S(V) 
for h — O.ber and for h = 5e F - In the former case for all 
the curves with 7 > 0.6er, the tunnelling rate due to the 
ferromagnetic lead is greater than the exchange energy, 
and no peak splitting appears. When 7 becomes less 
than the exchange energy, a splitting of the peak starts 
to appear. In Fig.^jb, due to a larger exchange energy, 
the zero-bias peak is shifted to finite bias eV = h and 
for small 7 the splitting of the two peaks is now 2e F re- 
flecting the shift to finite bias of the induced gap in the 
density of states (the gap is now that for one spin-species 
only). 

To compare with numerical results, we repeat the cal- 
culation used to obtain Fig]| but this time there is a finite 
exchange energy in the normal lead and central region. 
The results are shown in FigJ| where h = 0.025. We 
again see that we have good qualitative agreement with 
the whole curve being shifted to finite bias. 

Finally we show that this agreement can be made al- 
most exact. In Fig. [}] we show results for the case 7 ex- 
In this case equation (|^) can be approximated by (for 
A > e with r = 0) 
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We see that the sub-gap conductance is just a Loren- 
zian: at T = h = for A 3> e equation (^) is simply 
S(V) = S , (0)7 2 /(( e y) 2 +7 2 ). In Fig.|a we show the sub- 
gap conductance obtained numerically with tp'F = 1-0 
and tps = 0.2 averaged over 500 realisations. Each curve 
corresponds to various values of h = 0.00, 0.01, . . . 0.19. 
The value of the gap was A = 0.17. Each curve has 
been translated vertically by an amount 0.15 for clarity. 
The curve for h = can be fitted well by a Lorentzian 
plus a constant C to take into account the fact that in 
the numerical result we cannot neglect completely the 
other terms in equation (^J) . From this we obtain a value 
of 7 w 0.015 and 5(0) = 0.23 and the constant was 
C = 0.075. Putting 5(0) and 7 into the above formula 
at T = for the same values of h and A used in the 
numerics gives us Fig.^b. 

In conclusion, we have studied the sub-gap conduc- 
tance of a F'/F/S structure for varying ratio of the inter- 
face resistances. The presence of an exchange field yields: 



i) a shift to finite bias of the zero-bias peak when the F/S 
resistance dominates; ii) a splitting of the finite-bias peak 
when the two resistances are comparable. 
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FIG. 1. The system under consideration. 
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FIG. 2. The normalized conductance S = (dI/dV)/G n 
as obtained from quasiclassical theory, for varying 7, for 
j/er = 0.05 (lower most curve at large V), 0.25, 0.45, 0.65, 
1.0 and 2.5 (the uppermost curve at large V). The applied 
voltage is in units of the er- Here the exchange energy h = 0. 
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FIG. 3. The conductance G as obtained from numerical 
simulations, for varying 7, from 7 > er uppermost curve to 
7 < £t lowermost curve. Here the exchange energy h = 0. 
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FIG. 4. The same as in figure 2 for finite exchange energy, 
(a) h = 0.5et and (b) h — 5er- 
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FIG. 5. The same as in figure 3 for finite exchange energy, 
h = 0.025. 
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FIG. 6. (a) The conductance G vs voltage as obtained from 
numerical simulations, for varying exchange energy, h, from 
h — (lowermost curve) to h = 0.19 (uppermost curve). 
These results are in the regime 7 > et. Each subsequent 
curve was translated vertically by an amount 0.15. A value 
of A = 0.17 was used, (b) The normalised conductance S 
vs voltage as obtained from quasiclassical theory, for varying 
exchange energy, h, from h — (lowermost curve) to h = 0.19 
(uppermost curve). Each curve is shifted by 0.15 vertically. 
Values for the zero voltage conductance and gamma were ob- 
tained by fitting to the numerical result for ft = and the 
same value of A was used. 



5 



